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In previous papers! the authors have presented measurements of the 
critical absorption wave-lengths belonging to the L series of nine of the 
heavy elements (27 wave-lengths in all). We have also discussed the 
bearing of these data on theoretical matters, and their relation to other 
available measurements of X-ray spectra. Since these papers were 
published, many accurate measurements have been made of the emission 
spectra in both the L and M series for the heavy elements. 

Coster,? working in Siegbahn’s laboratory, has accurately measured 
the L series emission lines for eleven heavy elements. He has extended 
the three known M absorption limits to five for uranium and thorium, 
and has measured the first three M limits for bismuth. 

De Broglie and Dauvillier* have studied in detail the L series emission 
spectra for six heavy elements. 

With these new, accurate data, Smekal,‘ Coster,? Dauvillier® and 
Wentzel® have constructed atomic systems to explain the source of each 
known emission line. The systems are based on the Bohr-Sommerfeld 
theory of the origin of characteristic X-radiation. In developing and 
discussing their systems, these investigators have used the values of the 
L absorption limits measured by us. Certain criticisms have been made 
of a few of our values. It has therefore seemed desirable to measure 
some of the limits again. 

Coster® has suggested a possible source of error in our values for the 
absorption limits of tungsten on the ground that we used an X-ray tube 
equipped with a tungsten anticathode, and that therefore we would get 
an effect not due to absorption alone but also to the neighboring emission 
lines. We wish to call attention to the fact (stated on page 534 of our 
paper in the Physical Review, Dec. 1920) that we took’ the precaution 
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of using an X-ray tube with a molybdenum target to avoid exactly this 
difficulty. If there had been an appreciable amount of tungsten on the 
anticathode due to sputtering, we would have detected the presence of 
the strong emission lines of the L series. In a more recent paper Coster® 
states that he has measured the L; and L, absorption limits of tungsten and 
obtained essentially the same wave-lengths as those given by us. 

Table I contains our new values for the L, and L; absorption limits of 
Pt, Au and Bi. These measurements were taken with almost the same 
arrangement of apparatus as before. The X-ray tube was provided with 
an arm which reached nearly to the first slit of the spectrometer, thus re- 
ducing the absorption of X-rays by air. At the end of this arm and also 
on the front of the ionization chamber thin mica windows were employed 
which still further reduced the absorption of X-radiation. Finally the 
sensitivity of our electrometer had been increased many times. We 
were thus supplied with a much more sensitive detector of X-ray spectra 
than in our earlier measurements. This enabled us to use narrower 
slits and to obtain at the same time greater drops in our ionization cur- 
rents for the absorption limits. 

TABLE 1 


CrrricaL ABSORPTION. L SERIES oF X-Rays 
GRATING SPACE FOR CALCITE, 2d = (6.056 + 0.004) X 1078 Cm. 


CHEMICAL Le Ls 

ELEMENT A X 108 cm. */v@ XX 108 cm. */v@ 
Pt .9318 = 3 0.9779 .8918 = 9 1.022 
Au .9008 += 4 1.0116 .8610 + 9 1.058 
Bi .7871 + 4 1.1578 .7562 += 9 1.205 


In the first column for each limit is given the wave-length in Angstrom 
units; in the second, the wave number divided by the Rydberg constant. 

Our new value for L, of Au differs from our previous one by 0.17 per 
cent. The error in our previous work was due to the proximity of the 
bromine ionization limit, which has a wave-length = 0.91796 X 10-* cm. 
In our new measurements it was unnecessary to have ethyl bromide in 
the ionization chamber. 

For the L; limits to Pt and Bi the new values differ from the old ones 
by a considerable amount in excess of our then estimated error. These 
errors were undoubtedly inherent in the faintness of the effect which we 
were endeavoring to measure. 

The other three values are substantially the same as we gave before. 
We do not find the error in L; for Au suggested by Dauvillier. 

At the time of our previous investigation, the available evidence seemed 
to indicate that the 8; and 2 emission lines belonged, respectively, to the 
L, and L, groups—i.e., that 8; was emitted when an electron fell into the 
L, orbit, and 2 when an electron fell into the L, orbit. As we definitely 
proved that the L; and L, absorption limits for tungsten were of lower 
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frequency than §; and 2 (a fact since verified by Coster), the idea of 
absorption limits being “limiting frequencies’ for groups of emission lines 
was confronted with an anomaly. We proposed an explanation of. this 
apparent anomaly that was in accord with the Bohr theory. Since then, 
however, Hoyt® and Dauvillier® have pointed out that the W emission 
line at \1.212 should be taken as #8; instead of the one at 1.203 as before. 
This suggestion has been confirmed by the more accurate measurements 
of the corresponding line for the other heavy elements. When the square 
root of the frequency is plotted against atomic number for the. f; lines, 
the points all lie on a smooth curve. Furthermore, in place of 2, the work 
of Dauvillier and Coster indicates that a new y line, ye, at + = 1.072, 
is the result of an electron falling into the L, orbit, and therefore belongs 
to the L, group. Both the old {; line for W and the 72 lines have been 
assigned to the L; group. 

Consequently, we now take 8;, ys and +, to be the shortest emission 
wave-lengths in the L;, Ll, and L; groups, respectively. 

Coster has compared his new values for the wave-lengths of each of these 
emission lines with our values of the corresponding absorption limits, 
and has concluded that the absorption limits are in general somewhat 
shorter in wave-length than the emission lines associated with them. 
Table II gives this comparison using our new values for L, and L3. 


TABLE II 
ELEMENT In 6 Ls 1 
Pt .9318 = 3 .9317 .8918 = 9 .8950 
Au .9008 = 4 .90125 .8610 = 9 .8663 
Bi .7871 + 4 .7874 .7562 = 9 :7610 


Within experimental error the wave-lengths of the absorption limits 
are shorter than the wave-lengths of the shortest emission lines associated 
with them. 

Coster? now places the number of L doublets at four, n —1, B2—a2, ys— Bs, 
"e— Bs. ‘The y1—82 doublet has been eliminated which is in accord with 
our suggestion that the electrons responsible for these two lines do not 
appear to start from the same orbit but different ones. He has also found 
and measured the two M absorption limits predicted by us for U and Th 
in order to explain the existence of the n—/ doublet. His values are in 
close agreement with our predicted ones. 

According to Smekal,* the My line is due to the passage of an electron 
from an N; orbit to an M; orbit. If this is true, the following relation 
between frequencies should hold: 

(Ls; — Li) — My = LA — Lp 
Coster has tested this relation using his new values for 62 and f;, Sten- 
strom’s value for My and-our values of L; and L;. He finds it confirmed 
by the available experimental data except in the cases of Pt and Bi. Our 
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new values bring Pt into line. For Bi, however, our correction is in the 
right sense but far too large. We cannot at present explain this dis- 
crepancy. Table III gives our new values for this relation. 


TABLE III 
ELEMENT Le-Ih My (L3- Ii) -M B3:- 8 
Pt .1706 .17158 — .0010 .0000 
Bi .2163 .20144 + .0149 + .0176 


The lines 8; and y, have been assigned by Coster to electrons which start 
from different outer orbits and fall, respectively, into L, and L; orbits. 
Assuming the , electron to start from a higher orbit than 8;, we should 
then expect y; — (L; — Li) to be greater than 8;. Coster shows this 
to be the case except for Pt and Bi. Our new measurements bring both 
Pt and Bi into line with the other elements. Table IV contains values 
for these quantities in the cases of Pt and Bi. 


TABLE IV 
ELEMENT Ma— (L3- L1) Br A 
Pt 847.87 844.9 +2.9 
Bi 980 .6 (977 .2) +3.4 


The value of 8; for Bi was obtained, as suggested by Coster, by inter- 
polation based upon the constant difference in wave-length exhibited by 
the elements for B. — 8:, the wave-length for 6. being known. The 
differences in the fourth column are in the same sense and of the same 
magnitude as those found by Coster for Ta, W, Au, Tl and Pb. More- 
over, if we estimate a value for the wave-length of y, for U based upon 
our value of L; for the same element, we obtain a difference greater than 
those above and also in the same sense. It is therefore quite well es- 
tablished that the 8; and +, electrons start from orbits which maintain 
approximately the same difference in energy level as we go from Ta to U. 

The authors in 1920 compared the frequencies of certain L emission 
lines with the differences in frequency between corresponding L and M 
absorption limits.‘ Kossel first pointed out that these quantities should 
be equal on Bohr’s theory. At that time we found each L emission fre- 
quency to be equal to a corresponding difference in absorption frequencies 
within the limits of experimental error. Coster?* has recently extended 
this comparison with his new and more accurate measurements to in- 














TABLE V 
ELEMENT 1 Li--Ms n I:—M; a Li- Mi a Ii— Mz 
U 855 .84 855.1 1134.95 1131.8 1003.23 1002.17 990.37 989.21 
Th 819.19 818.8 955.78 954.80 944.08 943.15 
Bi 798.54 797.43 790.20 789.35 
a In—M:2 Bs L:— Ns Bs Ls— Mi 


U 1269.08 1265.87 1286.29 1285.78 1222.53 1220.86 
Th 1194.94 1193.15 1211.67 1209.77 1155.00 1151.86 
Bi 959.93 958 .16 973 .85 971.2 ——- ae 
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clude two more M absorption frequencies for U and Th, and the first three 
M limits for Bi. ‘Table V contains his values with the exception of L; — 
M; for Bi which we have corrected to accord with our new value of L; 
for that element. 

In each case the emission frequency exceeds the difference in absorp- 
tion frequencies, but by an amount that is less than 0.3 of a per cent. 

The K Series of Molybdenum.—The authors have also measured the 
emission and absorption spectra for the K series of molybdenum. ‘The 
apparatus and method used are described in previous papers.’ A Coolidge 
tube equipped with a molybdenum target and an extra arm at the end of 
which happened to be a thin glass window served as the source of ra- 
diation. * 

Table VI contains the results of these measurements. 

TABLE VI 


K SERIES OF MOLYBDENUM 
GRATING SPACE OF CALCITE = 2d = 6.056 (10) —® Cm. 


LIME NO. GLANCING ANGLE WAVE-LENGTH X (10)* cm. ORDER 
a2 3 13 °-36’-08” 0.71213 + 8 2d 
ay 3 13 °-31'-07" 0.70783 + 7 2d 
Bp’ 3 12°-02'-50" 0.6320 += 2 2d 
B 6 12°-01'-50”" 0.63114 + 7 2d 
BB’ 1 5°-59’-05" 0.6314 + 3 Ist 
7 3 11°-48’-32” 0.6197 += 2 2d 
¥ 1 5°-52'-30” 0.6199 = 3 Ist 
K 1 5°-51’-40” 0.6184 += 2 Ist 


In the second column the number is the number of actual measurements 
made. : 

In the second order, the 6 line was markedly asymmetrical and we were 
able to estimate the position of its weaker component, 6’, predicted by 
Sommerfeld. The wave-length difference between these two lines is of 
the same order of magnitude as that found in the case of other elements. 
The wave-lengths of 8 and + found in the first order are definitely longer 
than those found in the second order. This is chiefly due to the complex 
structure of both of these lines, although several investigators including 
the authors have found a similar but smaller effect in the case of emission 
lines supposed to be. single. 

In a recent paper Overn!® has also given wave-lengths for these same 
emission lines. His values, however, do not agree with ours even when 
corrected for the difference between the rock salt grating space used by 
him and the calcite grating space used by us. The discrepancy between 
his values and ours seems to increase as we go farther from the ‘white 
ine’ used as a reference line by him. 

As we determined these wave-lengths by the ionization method, we are 
also able to give the ratio of the intensities of a; and a. From three de- 











90 PHYSICS: CLARK AND DUANE Proc. N. A. S. 


tetihinations in the second order we obtain for this ratio 1.97 as compared 
with: 2.00 predicted by Bohr. 

The'ratio of the intensities of 8 to y in the first order is 6.3 (one deter- 
mination) and in the second order is 5.46 (two determinations). 


1 Duane and Patterson, Proc. Nat. Acad. Sci., Washington, Sept. 1920; Phys. Rev., 
Dec.‘ 1920, p. 526. ’ 

2 Coster, Z:.Phys., 5, 1921 (139); also 6, 1921 (185). 

3. De Broglie and Dauvillier, Paris, C. R., 173, p. 187. Dauvillier, Ibid., 172, also 
173, p. 647, also other papers in this volume. 

4 Smekal, Z. Phys., 5, 1921 (91, 121). Also 4, 1921. 

5 Wentzel, Fbid., 6, 1921 (2). 

6 Coster, Ibid., 5, 1921 (147). Also Physic Rev., Sept. 1921. 

7 Duane and Patterson, Physic Rev., 16, 1920 (534). 

8 Coster, Ibid., Jan. 1922. 

® Hoyt, Proc. Nat. Acad. Sci., Washington, Nov. 1920. 

10 Overn, Physic. Rev., Nov. 1921, p. 350. 





A- NEW METHOD OF USING X-RAYS IN CRYSTAL ANALYSIS 
By GEorGE L. CLARK AND WILLIAM DUANE 


JEFFERSON PuysicaL LABORATORY, HARVARD UNIVERSITY 
Communicated April 4, 1922 


The procedure described in this paper differs in some fundamental de- 
tails from the elegant methods of analyzing crystals by means of X-rays 
developed by Laue, the Braggs and others. 

The substance to be examined may be in the form of powder or of whole 
crystals. In the latter case the crystal need not be large. One having a 
volume of only .25 cubic mm. suffices, although slightly more accurate 
data may be obtained, if it has a somewhat larger volume. 

In our method we employ X-rays belonging to the continuous spectrum, 
and therefore containing rays of a great variety of wave-lengths, and not 
the line spectra of definite wave-lengths, as in most other methods. This 
procedure has certain advantages. An ordinary X-ray tube with a tung- 
sten target may be used, and high voltages may be applied to it. More 
intense X-radiation can be obtained from tungsten than from such chem- 
ical elements as rhodium and molybdenum. Further in some cases it is 
desirable to employ X-rays of shorter wave-lengths than those in the K 
series of rhodium and molybdenum, for these shorter rays are more pene- 
trating than the longer ones. As a rule we use the spectrum between .12 
and .80 Angstrém. 

In making an analysis of a small crystal we place it at the center of an 
X-ray spectrometer, provided with an ionization chamber and a quadrant 
electrometer. The slits of the instrument should be fairly narrow, and 
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the best results can be obtained by turning the X-ray tube so that the 
rays that pass through the instrument leave the target almost grazing its 
surface. This gives a narrow beam of rays. We set the ionization cham- 
ber so that the line from the axis of rotation to the chamber’s slit makes 
an angle, 20, of convenient size (usually about 10°) with the direct beam 
of rays. Keeping the ionization chamber always in the same position we 
turn the crystal around by successive small steps and measure the ioniza- 
tion current at each step. When a set of parallel planes in the crystal 
lie at an angle @ to the incident beam of rays the planes pick out of the 
beam X-rays of certain definite wave-lengths, , and reflect them into 
the chamber, and we get a marked increase in the ionization current. 
Evidently the angles between the positions of the crystal that corre- 
spond to maximum ionization currents are the angles between the crystal 
planes. 





e— 


FIGURE 1 


The curve in figure 1 represents the ionization current as a function 
of the position of the crystal in an experiment with potassium iodide. 
The salt came from Professor Baxter’s laboratory, and was of exceptional 
purity. The peaks on the curve are sharply marked, and their angular 
positions can be estimated to within less than 0°.01. The angles between 
several pairs of peaks are 44°.96, 18°.37 and 26°.52, which correspond very 
closely with the angles between some of the principal planes in a cubic 
lattice, namely, 45°, 18°.43 and 26°.57. This together with a similar ex- 
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periment performed with the crystal turned through 90° indicates that KI 
has a cubic lattice, which agrees with previous work by Davey and Wyckoff. 
The accuracy with which the experimentally determined angles agree 
with those given by solid geometry amounts in reality to a test of the 
precision with which the instrument has been set up, i.e., so that an axis 
of the crystal lies parallel to the axis of rotation, and so that the center 
of the incident beam of rays and the effective reflecting plane of the crystal 
lie on the axis of rotation. 

One peak of the curve, marked X, does not correspond to any of the 
chief planes of the crystal, i.e., those with small indices. Its position 
shows that the corresponding planes make an angle of 17°.84 with the 100 
planes. Possibly the marked reflection at this angle may be due to the 
characteristic line spectrum of iodine. Further evidence on this point 
will be given in a later paragraph. 

The wave-lengths, \, of the X-rays reflected at the angle 6, corre- 
sponding to any one of the peaks, are given by the fundamental law of 
X-ray reflection, which may be expressed by the equation 

nr = 2dsin 8, (1) 
where 7 is a whole number, and d is the distance between successive planes. 
We cannot use this equation alone to determine d, for we do not yet know 
the value of \. The peaks on the curve in figure 1 give us the angles be- 
tween various planes only. 

In order to determine the distance, d, we make use of the fact! that, 
in order to produce an X-ray of wave-length y, the voltage, V, applied 
to the X-ray tube must have a value at least as large as that given by the 
equation 

Ve = x? (2) 
where h is Planck’s radiation constant, e the charge of an electron (e) and 
c the velocity of light. This equation may be written in the form 


j 
Vi = — = 12,354, (2') 
@ 


where the numerical value’ of the universal constant, hc/e, corresponds to 
V expressed in volts and y in Angstréms. 

In an actual experiment to measure d we set the crystal so that the X- 
rays reflected by the planes we desire to investigate enter the ionization 
chamber (i.e., we set on one of the peaks in the curve of figure 1), and re- 
duce the voltage until the ionization current sensibly vanishes. The 
voltage at which it vanishes is sharply marked, especially if the voltage 
is a constant one, from a storage battery, for instance. It is a good plan 
to increase the current through the tube at the same time as the voltage 
decreases in such a way as to keep the product of the two, that is, the 
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power delivered to the tube, constant. Evidently the ionization current 
will not completely vanish until the voltage becomes so low that it cannot 
produce the longest X-rays reflected at the angle @ according to equation 
(1). The longest waves correspond to n = 1. 

A small correction has to be added to the angle @ due to the fact that 
the source of rays and the slit of the spectrometer have finite breadths. 
The rays in the incident beam are not all exactly parallel to each other 
and the correction equals half the angle between the rays that 
deviate most from the central ray on its two sides. This angle is the 
breadth of the source plus the breadth of the slit divided by the distance 
between them. ‘The correction is usually small. In our experiments it 
amounts to about 3’ of arc, which is less than one per cent of the corrected 
value of @ = 5°.51 used in the first illustrative case described below. 





FIGURE 2 


The curves in figure 2 represent the ionization currents as functions of 
the voltages for three of the peaks of figure 1. The voltages at which the 
ionization currents vanish for the peaks marked 100 and 110, respectively, 
are 18,600 and 26,310 volts. 

Eliminating \ between equations (1) and (2’) we get, with » = 1, 

_ ) 
2Vsin 6 
Substituting in this @ = 5°.51 and the values of V, we find for d the 
3.535 and 2.495 Angstroms, respectively. The ratio of the first to the 
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second is 1.417, which differs little from V2 = 1.414, the value it should 
have in a cubic crystal. Taking 3.532 as the average length of the edge 
of a cube, and assuming one atom to éach cube the calculated density of 
KI comes out 3.111, in close agreement with the density 3.114 measured 
by Baxter. 


Equation (3) contains a universal constant— = 12,354 and two 
e 


quantities, a voltage, V, and an angle, 6, both of which we measure. . 
The values of d, therefore, determined by this method are independent 
of previous measurements of the distance between the planes in any 
crystal. This advantage is, perhaps, more theoretical than practical at 
present, for measurements with X-rays reflected from calcite probably 
furnish us the most accurate estimates we have for the universal constant. 





Kilovelts —> 
FIGURE 3 


The agreement between the distances, d, determined by this method with 
those calculated from the density, atomic weights, etc., helps to confirm 
the general correctness of the laws of X-rays as at present understood. 
If the specimen to be examined has the form of powder, we place the 
powder at the center of the spectrometer, fix the position of the ionization 
chamber at an angle 26 and reduce the voltage by steps making measure- 
ments of the ionization current as before. The ionization current vanishes 
when the voltage is no longer able to produce the longest X-rays that are 
reflected at the angle @ by the powder. This maximum wave-length 
corresponds to m = 1 in equation (1) and to the greatest distance d be- 
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tween successive planes in the crystal that are alike. The curve in figure 
3 represents the ionization current as a function of the voltage in an 
experiment with potassium iodide powder. In addition to the point 
where the ionization current. vanishes there are several well-marked 
places where the curvature changes abruptly. The voltages at these 
breaks in the curve are the voltages at which the rays reflected from other 
crystal planes cease to be produced. These voltages together with the 
angle 6, substituted in equation (3), give the corresponding distances be- 
tween the planes. The three most sharply determined voltages are 
21,740, 31,010 and 37,600 volts, respectively. The ratios of these numbers 
to each other differ from V2 and V3 by small fractions of one per cent 
only. They correspond, therefore, to the 100, 110 and 111 planes, re- 
spectively, of a cubic lattice. The angle @ equals 4°.51 in this experiment, 
and the average value of d, calculated from equation (3), is 3.534, in close 
agreement with that obtained in the above-described measurement with 
a small crystal of potassium iodide. 

The curve marked X in figure 2 corresponds to the peak marked X in 
figure 1. There appears to be a sharply indicated point of inflection in 
the former at a voltage of 30,000 volts. Substituting this in equation 
(2") we find that the corresponding wave-length is .412. This wave- 
length is a good mean value of the wave-lengths of a and £ lines in the K 
series of iodine, namely, .437 and .388. 30,000 volts do not suffice to 
produce the K series of iodine. If, however, one sets a spectrometer on 
a characteristic line and reduces the voltage, the ionization current does 
not vanish when the voltage ceases to produce the characteristic X-rays. 
It vanishes only when the voltage becomes insufficient to produce an X- 
ray in the continuous spectrum of equal wave-length to that of the char- 
acteristic line. We have tested the existence of the peak X by turning 
the crystal as in the experiments represented in figure 1 with voltages 
in the neighborhood of 33,030 volts, which is the critical voltage for the 
production of the K series of iodine. We find no evidence of a peak at a 
voltage slightly below 33,000. The peak is just observable, however, 
if the voltage equals 33,500 volts. These points strongly favor the idea 
that the characteristic K radiation of iodine produces the peak X. The 
distance between the planes making an angle of 17°.84 with the 100 planes, 
however, does not equal that calculated from equations (1) and (3). The 
former is 1.088 and the latter is 2.144, i.e., about twice the former. Per- 
haps this may be due to the fact that only half the atoms are iodine. 
If the characteristic radiation does produce the peak, and, if we adopt 
the interference of scattered rays as the explanation of X-ray reflection, 
we must assume that, on the average at least, definite time relations must 
exist between the passage of the primary X-ray across an atom and the 
emission of@ secondary characteristic ray. If we adopt the quantum point 
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of view, we must suppose that the direction in which the characteristic 
ray travels depends upon neighboring atoms as well as upon the angle of 
incidence of the primary beam. 

The relative reflecting power of different sheen may be determined by 
reflecting from them the same portion of the continuous spectrum, or differ- 
ent portions having known relative intensities. The use of the con- 
tinuous spectrum for this purpose instead of the line spectrum has the 
advantage of flexibility. One can employ short, penetrating X-rays, and 
with these investigate dense crystals of high absorbing power. 


1 Phys. Rev., Ithaca, Aug. 1915, p. 166. 
2 J. Optical Soc., May 1921, p. 386. 





ACOUSTIC TOPOGRAPHY VARYING WITH THE POSITION OF 
THE ORGAN PIPE* 


By Cari Barus 
DEPARTMENT OF PHYSICS, BROWN UNIVERSITY 
Communicated March 20, 1922 


1. Introductory.—The coérdinate axes have the x-y plane on the large 
table, the origin near one corner. The survey is made between walls 
at y = 174 cm. and y = —130 cm.; and the pipe is displaced both in 
the y and z (upward) direction. In the graphs, the letters P, R, T, T’, 
used either as subscripts or directly, will refer to the pipe, the reso- 
nator, the edges of the table. Throughout s is the fringe deflection of 
the U-tube interferometer, proportional to the acoustic pressure or to 
the nodal density evoked at a given place (x = 0, y, 2), by the sounding 
pipe P. 

The latter was a closed brass f’’ pipe, blown by the pipe blower else- 
where described, rigidly attached. In graphs 1-3, this is sketched 
in under P and an arrow shows the direction of the wind current. The 
closed cylindrical pinhole resonator R, lying with its axis parallel to 4, 
in the yx plane (table) and moving end-on along y, is also suggestively 
indicated, showing the resonator azimuth for which the usually paired 
graphs apply. In the earlier work (figures 1-3) the middle of the resona- 
tor was used for the location of points in y. An inversion of the reso- 
nator thus gave the two curves for each case and they were obtained con- 
secutively. If the mouth of R is used to determine the codrdinate y 
as in figures 4, 6 to 10, but a single curve is obtained, for all azimuths of 
R. ‘The same is true for a resonator placed with its mouth at the mean 
Yr of a given ordinate, s, in case of figures 1-3. 

2. Pipe Elevated above the Origin on the Table—In figures 1-4 the or- 
_ gan pipe, with its mouth vertically above the origin, was ragged succes- 
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sively in steps of 20 cm. from z = 20 cm. to z = 80. For each pipe level 
the distribution of nodal intensity along the y axis was explored by mov- 
ing R along it (with observations taken every Ayer = 10 cm. apart). In 
figure 1, z, = 20 cm., there is probably some distortion at the origin, 
owing to the wind from the near pipe, otherwise the first crest would 
be in negative y. From the glancing incidences, this graph has less 
detail than figure 2 (z, = 40 cm.), which preserves the same general 
character but is lower in s and shrunk in the y direction. In figure 3, 
z, = 60 cm., the topography is totally changed. There is now a trough 
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near the origin and (on compounding) but one marked crest. Figure 
4, 2» = 80 cm., is redeveloping a crest at the origin, but otherwise pre- 
serves the former outlines. In this figure and the following, the mouth 
of the resonator was used to define yz so that the two curves are practically 
coincident. 

If the highest crests correspond to like path differences and there- 
fore belong together, their position yp, 2, in figure 5 is seen to progress 
roughly along a straight line y = 1.252 — 40, the slope dy/dz corre- 
sponding to an angle of incidence of a = 51°. The minima determine 
parallel lines, the vertical distance apart of consecutive lines being here 
Az, = 24 cm. or about half a wave-length. ; 

Furthermore, the graphs, as a whole, are alternately high and low in 
intensity, s. If the intensity s,, of the highest crest be expressed in terms 
of the pipe elevation z,, the resulting graph, figure 2’ (scale '/2), presents 
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the same strongly marked harmonic character of the reflection graphs 
of the preceding paper. Cf. § 3. 

3. Pipe Displaced in a Given Level along y.—Figures 6-10 are constructed 
on the same plan as the preceding, the ordinates s, denoting fringe de- 
flections or the nodal intensity, encountered by the mouth of the res- 
onator at yr. If the ordinates be also taken as 2,, the pipe position, 
P, etc., may also be indicated on the graphs. TJ, 7’ show the edges of 
the table. In figure 6, where the pipe is 40 cm. above the origin, the graph 
is highly detailed with crests mostly 30 cm. apart. The hyperbolas 
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for reflections from the wall at y = 174 cm. and beyond T’ have been 
sketched in at 5d and 6A of path difference, together with the corre- 
sponding acoustic rays. So also those for reflection from the wall at 
y = — 130 cm. (beyond 7). The last agrees badly with the crest observed 
at about —10 cm. The alluring possibility of an intersection of the 
positive and negative hyperbolas resulting in very high crests does not 
work out. In figure 7 also the positive hyperbola 5) fairly corresponds 
to the observed crest, but not the negative one. In figures 8, 9, 10 the 
positive hyperbolas 3\, 4A, 5A, also, fail of coincidence with the observed 
crests. These in fact lie nearly between them, or half a wave-length of 
path difference beyond. Thus it seems improbable that the maxima and 
minima can be explained as resulting from the reflection at these walls, 
while all other walls and the ceiling are far more distant. 

In figures 6-10 one is often in doubt as to which of the crests belong 
together. Crests seem rather to separate and recombine. Taking the 
highest in the successive curves, figure'11 shows that the relation of pipe 
and resonator positions for maxima is here not simple as in figure 5. In 
fact the pipe approaches the high crest and finally overtakes it (figure 10). 
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There is no symmetry of graph at the middle of the table (y = 33) and 
the two edges T, T’ show no correspondence. ‘The graphs may actually 
be traced into free air beyond the table, as indicated in the prolongation 
of figure 12. Finally if the intensity s,, of the supposedly chief maxima be 
laid off in terms of the pipe position an harmonic figure 12 of relatively 
constant amplitude appears in reasonable contrast with the curve of di- 
minishing amplitude in figure 2’, for the raised pipe. In both cases, the 
steps of 20 cm. between pipe positions are too large for sharp discrimina- 
tions; but it is noticeable that the distance between crests in figures 2’, 12 
and the earlier reflection figures is about Ay = 40 cm. ‘The prevailing 
Ay in the single graphs 2, 6, etc., is about 30 cm. In figure 12 the points 
beyond y, = 80cm. were obtained in similar investigations, here omitted. 

4. Summary.—The above work as a whole has shown that when the 
reflecting plate is relatively near, to the pipe and resonator, the position 
of maxima and minima may be predicted as a case of ordinary interfer- 
ence and the wave-length computed satisfactorily. The distribution of 
intensity among the crests and troughs remains harmonic and has not 
been foreseen. 


If the reflecting surfaces are relatively remote the positions of crests 
and troughs cannot, as a rule, be found by the same method satisfactorily. 
In certain instances, there seemed to be an approximate fit; but as a whole 
the attempt was unsuccessful and the distribution of nodal intensities 
equally puzzling. Acoustic topography is not symmetrical to the pipe. 
It appeared, however, that, when the organ pipe, definitely placed, is 
sounded, the occurrence of nodal surfaces of a fixed position in free air is 
demonstrable everywhere. Distributions between walls are consistent 
and differ from distributions between wall and door. It is not improbable, 
moreover, that such surfaces are regularly grouped and may for reasonable 
distances be approximately parallel. If then they are intersected obliquely, 
for instance by the plane of the table, and if a is the mean angle be- 
tween the latter and the nodal surface, the distance between crests on 
the table should be Ay = 'd/(2 sin a). Hence this intercept, Ay, may 
have any value depending upon the general shape of the room. If the 
prevailing value Ay = 30 cm. be taken, sin a = 24/30 or a = 53° 
roughly. This angle between surfaces is also the angle of incidence of the 
rays, and has been encountered more or less closely in other instances 
(fig. 5). The distribution of intensities would then also depend on the 
fixed pipe position in relation to the given room, regarded integrally as an 
enveloping reflector. Such an explanation is plausible and flexible, no 
doubt; but it has the great disadvantage that none of the results can be 
reproduced by computation. The conviction that some other explana- 
tion may be found is not removed and the answer probably lurks in such 
curves as figures 2’ and 12. This is particularly so, because the pin- 
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hole resonator is practically unresponsive to antinodes and to wave trains 
except within ranges from the pipe which are small as compared with the 
dimensions of the room as a whole. 


* Advance note from a Report to the Carnegie Institution of Washington, D. C. 





ON STEERING AN AUTOMOBILE AROUND A CORNER 
By A. G. WEBSTER 
DEPARTMENT OF Puysics, CLARK UNIVERSITY 
Read before the Academy April 25, 1921; received April 6, 1922 


The question of so arranging the steering gear of an automobile that 
the front wheels shall skid as little as possible owing to axes not inter- 
secting has been dealt with by A. L. Candy (American Mathematical 
Monthly, June 1920) and leads to a geometrical discussion of a high 
degree of complexity. The present somewhat trifling paper has no such 
ambitious purpose, but was suggested by a question asked on the witness- 
stand as to how much an automobile must turn out in order to pass another. 
It was necessary to make some assumption, and this led to the question 
of steering in general, and the principles that should govern it. We shall 
consider only the question of passing from a straight course to another at 
right angles therewith at constant velocity, and inquire how we must 
steer in order not to skid. 

The layman ignorant of geometry and dynamics will suppose that we 
describe a circle, but this is impossible, not only because it is impossible 
to put the helm over suddenly so as to change from zero to a finite curva- 
ture, but also because the sudden change in the centrifugal force from 
zero to a finite quantity would produce a wholly intolerable shock. In- 
asmuch as the question of transition railroad curves has received so much 
attention, it is thought that the treatment of this question for an auto- 
mobile may have some interest to engineers. We shall for the present 
neglect the fact that the rear wheels do not track with the forward ones, 
taking that up later as a correction. Thus we shall examine the curve 
described by a machine of zero wheel-base, like a wheelbarrow, measur- 
ing wheel or monocycle. The condition is evidently that the curvature 
shall begin at zero, increase to a maximum until the tangent has turned 
through an angle of 45 degrees, and then return by a symmetical process. 

The centrifugal force is equal to mv?/p, where m is the mass, p the radius 
of curvature, v the velocity, and if there is to be no skidding this must be 
less than mgy where uy is the coefficient of friction, consequently we must 
have at all times 


v?/p < mg. 
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This does not, however, determine the shape of the curve, to obtain which 
some assumption must be made. If we assume that the mean value of 
the centrifugal force or, what is the same thing, of the curvature (since v 
is constant) is to be a minimum we shall get no result, except a straight 
line, which would be a discontinuous solution, and inadmissible. We shall 
therefore assume that the mean square of the curvature is a minimum, 
that is 
5S nds = 5 f (d0/ds)*ds = 0, 
subject to the condition 
JS dx = fcoséds = constant, the length of chord. 


But this gives 

d*6/ds* + dsind = 0, 
where \ is some constant. But this is the equation of the elastica or 
curve assumed by a bent flexible rod, ruler, or spline. In fact, in order to 
obtain the equation of the elastica, we have to make the potential energy 
of bending a minimum, and as this, at each point, is proportional to the 
square of the curvature, the problem is mathematically the same. 

If we take the chord for the X-axis, we shall find the equation satis- 
fied if we take the curvature proportional to y, and introducing the elliptic 
functions of Jacobi, with c a constant of homogeneity, 

dx/ds = cos), dy/ds = siné, 


«x = d0/ds = —y/c?, dx/ds = d*0/ds? = — 1/c*dy/ds = — 1/c*siné, 
4 = acens/c = acnu, dy/ds = — a/csnudnu = — 2ksnudnu, 
whiere 


s/e= 4, a/c = Qk, 
from which 





dx/ds = V1 — 4ksn2udn’ = 2dn2u — , 
s/c s/c 
x = 2c fdn2udu — cf du = cl(2E(u)—u] = c[2E(k,e) —F(k,¢)], 
0 


where F and E are Legendre’s elliptic integrals of the first and second kind, 
respectively. The value of k is determined by the fact that at the point 
of inflexion, where 9 = —45°, we have «= K, the complete elliptic inte- 
gral, and since snK = 1 and dnK = V1—k, we have 2kvV1—k, = 
sin 45°. Putting k = sina we see that a = 22°.5. With this value for a 
taking values of g at intervals of 10° from 0° to 90°, values of x and y 
are calculated by means of Legendre’s table and the curve drawn, which 
is shown on the inside of figure 1. It fits very accurately a'spline made of 
celluloid, held fast at the vertex, and merely by contact at the point 
of inflexion. We find the ratio of the maximum ordinate to the whole 
chord or base is .2757 and that the radius of maximum curvature is the 
chord divided by 2.1429. 
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The question now arises, how nearly does the curve actually described 
correspond with the above description? We may take as the simplest 
assumption that when the driver gets ready to turn he puts the steering 
wheel over with a constant velocity until he is half-way around, and then 


ge 











/ 


FIG. 1 


turns it back with the same velocity. This is to be sure not exactly true, 
but as the wheel is dominated rather by friction than by inertia it will do 
for an approximation, and is besides about what observation shows. 
If the effect of this is to make the curvature of the path proportional ‘to 
the time or to the distance travelled since turning (which will be shown 
below to be nearly the case), we have 
x = d6/ds = as 

as the differential equation of the curve, where a is a constant. But a 
curve whose curvature is proportional to the length of the arc is the well- 
known spiral of Cornu, used in the theory of diffraction in optics. The 
curve may be constructed graphically from the differential equation, but 
more accurately by the use of tables of Fresnel’s integrals. If we put 
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x= fcosr/wdv, y = fsinr/ dv, 


we have 
dx = cosr/v'dv, dy = sinr/ dv, ds? = dx? + dy? = dy’, 
dy/dx = tan®? = tanr/.v*, 6 = r/? = r/2s? = as?/2. 


The integrals above have been tabulated by Gilbert, and are reproduced 
in Verdet’s Traité d'Optique Physique. From them the spiral has been 
plotted, and is shown in the outer curve in figure 1. It will be seen that it 
differs very little from the elastica. As a matter of fact the ratio of 
the maximum ordinate at right angles to the chord is .2861 compared 
with .2757 for the elastica. For the circle this ratio would be only .2071. 
The reason for the near agreement of the two curves, in one of which the 
curvature is proportional to the arc, and in the other to the ordinate, 
is of course the fact that for a considerable distance the curve is nearly 
straight. 


The simplest assumption that could have been made for a rough cal- 
culation would have been that of a parabola, y = ax’, but the ordinary 
parabola of order two would not do since its curvature does not vanish 
at the origin, and diminishes as we leave the origin. A cubical parabola 
would be the most obvious, but we shall do better by putting y = ax” 
and determining the exponent » so that the curve shall pass through the 
vertex of the elastica and be tangent to it. We find =3.458. The curve is 
drawn dotted on figure 1, and lies between the other two curves. “The ratio 
of the extreme distance from the chord to the length of the chord in the 
case n=3 would have been .25. We consequently see that any of the 
four curves, the elastica, spiral of Cornu, parabola of order 3.458, 
or cubical parabola, will give a very good approximation to the actual path. 

We have now to consider the effect of the finite length of the wheel- 
base, and this leads us to consider the question of steering in general. 
We shall first suppose for generality that both the front and rear wheels 
can be steered, as in a hook and ladder truck, but we shall still neglect 
the effect of the finite width of the track. The case contemplated is ex- 
actly realized in the bicycle. In figure 2, let / be the length of the wheel- 
base, let s, 0, a, p be the distance run, the angle made by the tangent with 
a fixed direction, the angle through which the wheel has been turned (that 
is, the angle made by the tangent to the path with the chassis, and the 
radius of curvature of the path, respectively, for the front wheel, and let 
the same letters with accents denote the corresponding quantities for the 
rear wheel. (When we come to consider the dynamics of the matter we 
shall need corresponding quantities for the center of mass, and shall de- 
note them by two accents.) Then since the whole machine turns about 
the instantaneous center which is the intersection of the axes of the two 
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wheels, we see on reference to figure 2, y being the angle made by the chassis 
with the fixed direction, 
pd@ = ds, p'dd’ = ds’, 0+ a=0' +a’ =y. 
We also have 
cosa/p = cosa'/p’ = sin(a—a’)/I, 
which equations give p and p’ as functions of a, a’. If then a and a’ 
are given as functions, say of s and s’, this makes seven equations between 











te] 
FIG. 2 


the variables a, a’, 0, 6’, s, s’, p, p’, so that we may find the differential 
equations of the path of either the front or rear wheel. In particular, if a 
and a’ are constant and equal, we shall have merely a movement of transla- 
tion, while if a’ = 0, a = constant we shall describe two concentric circles. 
If a’ = 0,a = 90°, the machine will describe a circle about the rear wheel. 
If, with these values of a, a’ the length / be made variable, and a given func- 
tion of s, the front wheei will describe a curve and the rear wheel its evo- 
lute. Thus a curve may be mechanically constructed from its intrinsic 
equation. ‘This principle is used in certain integrating machines, and the 
author has invented an integrator for drawing trajectories involving this 
principle, which will be described later if he succeeds in having it built. 
In the case of the automobile we have 

a’ = 0, p = lIctna, ds = Ictnadé, 
and if as before we assume the mode of steering to be a = as we may 
integrate and obtain the intrinsic equation of the curve, 

ld@ = tan(as)ds, la@ = log sec as. 

Developing the tangent in series we may see how much this differs from 
the Cornu spiral. The curve has been drawn graphically from the differ- 
ential equation, but since on a sixty-foot chord it differs but slightly from. 
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the spiral (which may be judged from the fact that for @ = 45°,a = 23°.2’; 
p = 22.35 compared with p = 22.28 for the spiral, and that the rear 
wheel’s track is less than two feet distant at the maximum) it is not here 
reproduced. 

It may be said that the above is not applicabie to the way in which an 
automobile actually skids, when it is well known that the two wheels do 
not break the frictional constraint together, but the rear wheel generally 
skids first. This may be true, but I think this will not invalidate anything 
that I have said. One would suppose that this question would be treated 
in any elementary book on statics, but as I have failed to find it in any 
treatise, I insert a treatment here. It will now be necessary to consider 
a third path, namely, that of the center of gravity of the machine, which 
we shall still assume to be on the line joining the wheels, at distances d and d’ 
therefrom, figure 3. Then we have d’ctna” = (d + d’)ctna. The three 





Oo 


FIG. 3 


paths have the same center of curvature. The constraints due to the 
friction can be only at right angles to the paths of the wheels, and we will 
call them F, F’. The kinetic reaction of the machine is described by the 
centrifugal force C, the tangential reaction R, and the moment S. It 
may be asked how the wheels can cause any reaction in the path, but it 
will be seen that there must be such when we consider that energy of 
translation being converted into energy of rotation and conversely will 
give rise to acceleration in the path, negative or positive, respectively. 
In fact we see at once that S must vanish with R, for three concurrent 
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forces do not give rise to any moment. Resolving parallel and perpendic- 
ular to the chassis, 
F’ + F cosa = Ccosa” + R sina”, 
F sine = C sina” — Rcosa’, 
and taking moments about the center of mass, 
F'd’ — Fdcosa = S, 
or making’ use of the relation between « and a”, 
F’ — Fsin(a—a”)/sin a” = S/d’. 
Solving’ these equations we obtain 
R = — S/d'sina’, 
F = sina”/sina (C + Scosa”/d’), 
F’' = sin(a — a”)/sina C + S/d’ (sin(a—a")cose”/sina —1), 


Now since C = mv?/p” and S = mk*d*/dt? = mk?(d?0/dt? + d?a” /dt?) 
and in the method of steering assumed this is mk*d?0"/di? = mkvd0/ds?, 
if v is constant, both terms are proportional to v?. Consequently the con- 
clusions above stated are justified. (It is to be remarked that if v is con- 
stant, the kinetic reaction above is to be replaced by the force of traction 
necessary to maintain the velocity constant.) The term in C is propor- 
tional to the curvature, and the term in S to the rate of increase of cur- 
vature, the latter is greatest at the beginning and end of the curve, the 
former at the middle. 

Finally I have made some experiments to determine the coefficient of 
friction » of rubber on slimy stone, with the conclusion that it is about one- 
fourth. Using the criterion v?< ygp we arrive at the conclusion that the 
safe speed for curves as above is somewhat less than ten miles per hour, 
which may account for the common traffic rule of ‘eight miles around the 
corner.” 





THE PROBABLE ERROR OF THE VITAL INDEX OF A POPU- 
LATION' 
By Joun RicE MINER 
DEPARTMENT OF BIOMETRY AND VITAL STATISTICS, JOHNS HOPKINS UNIVERSITY 
Communicated February 14, 1922 ; 
In comparing the birth-death ratios of. different communities we need 
their probable errors to determine the significance of observed differences. 


Let N = the population of a community, 
B = number of births in a year, 


D = number of deaths in a year. 100 B 
Then the birth-death ratio or vital index is defined as ———. 


D 
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The standard errors of sampling of births and deaths will be: 


=a SAS is D) 
63 = see a op = VY ———_—_ 


; 1— B/N 1—D/N 
hence w= 22, afin BIN rae 
= ona B 





and the probable error of the birth-death ratio will be? 


B.j1-—B/N 1—D/N 
oe 100 B = 67.449 5 V B + D (i) 
D 








since r = 0, being the correlation of errors of sampling of births and deaths. 
As B and D are small as compared with N we may write this approximately: 





B y! i 
“D 
This form has the advantage of not involving N, so that if the numbers 
of births and deaths in a community are known but not the population, 
the probable error of the birth-death ratio may still be found. 


For Baltimore in 1917 


N = 594,637 
B = 15,040 
D = 11,355. 


Therefore by (i) 





_15040 1 11355 




















1— _—_——_ 
a — 67.449 15040 V 594637 594637 
eee 11355 15040 11355 
747 
= 80.34 af 074T = 89.34 V.00006481 + .00008638 
+ See 
= 89.34 X .01229 = 1.10. 
By (ii) 
—— Va + —— = 80.84 V. 
PE (100 » = 67.449 Se eS +—=- = 89.34 V.00006649 + .00008807 


= 89.34 X .01243 = 1.11. 


In using these formulae it should be borne in mind that the expression 
for the standard deviation of an index, which has been used in their der- 
ivation, is obtained on the assumption that deviations are small compared 
with the mean values of the variables. This will be true if B and D are of 
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fair size, but the formulae should not be used for very small absolute 
numbers of births and deaths. 

1 Papers from the Department of Biometry and Vital Statistics, School of Hygiene 
and Public Health, Johns Hopkins University, No. 52. i 

2 Pearson, K. On a form of spurious correlation which may arise when indices are 
used in the measurement of organs. Proc. Roy. Soc., 60, 1897 (489). 


THE AVERAGE CHEMICAL COMPOSITION OF IGNEOUS ROCKS 
By FRANK W. CLARKE AND HENRY S. WASHINGTON 


U. S. GsoLocicaL SuRVEY* AND GEOPHYSICAL LABORATORY OF THE CARNEGIE 
INSTITUTION OF WASHINGTON 


Communicated March 22, 1922 


We have recently completed computations of the average chemical 
composition of the igneous rocks of various countries, of the continents, 
and of the earth’s crust. These are to be published as a Professional 
Paper of the U. S. Geological Survey, with discussions of the characters 
of the various averages, of the relative abundance, distribution, associa- 
tion, and evolution of the elements, and other matters. Inasmuch as it 
is probable that this paper will not appear for a considerable time, it is 
thought best to publish now the various average compositions arrived at, 
in order to put before geologists and others some interesting data as to the 
chemical composition of the earth. 

The data used by us were 5,159 chemical analyses of igneous rocks from 
all over the earth, and published between 1884 and 1913, inclusive. Analy- 
ses only of “‘superior” quality, and of fresh unaltered rocks, were con- 
sidered. These were taken from a collection of rock analyses made by 
one of us.' In computing the averages of the various countries and the 
continents the sum total of each constituent was divided by the total 
number of analyses of rocks from the area dealt with. In the case of 
the earth average, both in terms of oxides and of elements, the same 
method was adopted for the major constituents, SiO2, Al,O3, Fe2O;, FeO, 
MgO, CaO, Na,O, K2O and H,O (only water above 110° was considered), 
while the average for the minor constituents was obtained by taking the 
mean of the sum total of each divided by the number of analyses (which 
gives too low a result) and by the number of determinations (which is 
apt to give too high a result). 

In computing the average of the whole crust the relative proportions 
estimated by one of us? some years ago were adopted. The latest esti- 
mates of the masses and compositions of the hydrosphere and atmosphere 
have been used. These relative masses are as follows: Lithosphere 93%; 
hydrosphere 7%; atmosphere 0.03%. 
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AVERAGE IGNEOUS ROCKS, NORTHERN NORTH AMERICA 
1 2 3 5 
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. North America, including Greenland, Mexico, Central Aislin: and West Indies 
1709 analyses. 

. Greenland, including East and West Greenland and Ellesmere Land. 41 analyses. 

. Alaska. 24 analyses. 

. Canada, including Alaska, East Canada, and British Columbia. 188 analys 

. East Canada, including Ontario, Quebec, Maritime Provinces, and Newionaiiand. 
99 analyses. - 

. British Columbia, including Alberta. 60 analyses. 
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AVERAGE IGNEOUS ROCKS, MIDDLE AND SOUTH AMERICA 
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1. Mexico and Central America. 47 4. British Guiana. 45 analyses. 
yses. 5. Andes Mountains. 56 analyses. 
2. West Indies. 82 analyses. 6. Eastern Brazil (including Paraguay 
3. South America. 138 analyses. and eastern Argentina). 20 an- 
alyses. 
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PETROLOGY: CLARKE AND WASHINGTON  Préc.N.A.S. 


RARER ELEMENTS IN IGNEOUS ROCKS (Continued) 
1 
-000,0x 
-000,00x 
.000,00x 
-000,00x 


1. Clarke and Washington, 1922. 
2. Vogt, 1898. 


ELEMENTS IN THE LITHOSPHERE, HYDROSPHERE, AND ATMOSPHERE 
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100.000 . 100.000 100 .000 
1. Average composition. Ten-mile crust, hydrosphere, and atmosphere. 
2. Average composition. T'wenty-mile crust, hydrosphere, atmosphere. 
3. Average composition. ‘Ten-mile crust, igneous and sedimentary rocks. 
4. Average composition. Ten-mile crust. Igneous rocks. 
The lithosphere is assumed to be made up as follows: Igneous rocks 


95.00%; shale 4.00%; sandstone 0.75%; limestone 0.25%. 
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In the calculations involving the atmosphere and hydrosphere the thick- 
ness of the crust is assumed to be ten (or twenty) miles. 


* Published with the permission of the Director of the U. S. Geological Survey. 

1H. S. Washington, “Chemical Analyses of Igneous Rocks,” U. S. Geol. Survey, 
Prof. Paper 99 (1917). 

2 F. W. Clarke, ‘‘The Data of Geochemistry,” U. S. Geol. Survey, Bull. 695, 1920 (22 
and 33). 





DARK NEBULAE 
By Henry Norris RUSSELL 


RESEARCH ASSOCIATE, MouNT WILSON OBSERVATORY, CARNEGIE INSTITUTION OF 
WASHINGTON 


Communicated March 14, 1922 


1. It is now generally believed that many of the dark markings in the 
.Milky Way, and dark starless regions in the sky, are produced by the in- 
terposition of huge obscuring clouds between us and the remoter stars. 
A long list of such dark markings has been given by Barnard,' who has 
done more than any one else to point out their importance and probable 


nature. In some cases, as in the Pleiades, Orion and Ophiuchus, these 
“regions of obscuration’”’ merge into faintly luminous nebulosity in the 
vicinity of certain stars, in such a way that there can be no doubt that 
they lie near these stars in space. 

It thus appears that the obscuring masses or dark nebulae in Ophiuchus 
and Scorpius are at a distance of from 100 to 150 parsecs, those in Taurus 
at probably about the same distance, and those in Orion some 200 parsecs 
from us, while the dimensions of the individual clouds are themselves 
measured in parsecs. 

The occurrence of these three great regions of obscuration within a 
distance which is so small compared with that of the galactic clouds in- 
dicates that such objects are probably of great cosmical temperature. 

2. These dark nebulae usually appear to be quite opaque. In some 
cases the stars can be seen faintly through them, apparently without much 
change in color; but in certain instances? stars imbedded in dense luminous 
nebulosity are abnormally red. 

Of the various forms in which matter may be distributed in space, by 
far the most efficient in producing obscuration is fine dust, since this has 
the greatest superficial area per unit of mass. In a cloud composed of 
spherical particles of radius r-and density p, distributed at random so that 
the average quantity of matter per unit volume is d, the extinction of a beam 
of light in passing through this cloud will be e stellar magnitudes per unit 
of distance, where e¢ = 0.814 gd/pr. 
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The numerical factor is independent of the physical units which are 
employed. The factor q is introduced to take account of the complica- 
tions which occur when the size of the particle becomes comparable with 
the wave-lengths of light.* For particles more than two or three wave- 
lengths in diameter g is sensibly equal to unity. For smaller particles 
it increases and is a maximum, 2.56, when the circumference of the parti- 
cle is 1.12 times the wave-length. It then rapidly diminishes and becomes 
nearly equal to 14/3 X (2mr/d)‘ for particles of less than half this diameter. 
The ratio g/r is a maximum, 2.42, when the circumference equals the wave- 
length. 

For clouds of the same mean density d the opacity reaches a sharp max- 
imum when the particles are of this size. At the same time the absorption 
changes from the non-selective type to the selective type, varying as 
\~*. For visual light the maximum opacity occurs when the radius is 
0.086 p. A cloud of particles of this size, and of the density of rock (2.7), 
will exert an absorption of one magnitude if it contains only 1/86 of a milli- 
gram of matter per square centimeter of cross-section, regardless of its 
thickness. If the particles are of half this size, or smaller, the selective 
absorption is almost as complete as for a gas, but may be nearly 100 million 
times as great. 

Obscuration of light in space, therefore, whether general or selective 
with respect to wave-length, will be produced mainly by dust particles 
a few millionths of an inch in diameter, unless such particles form a negli- 
gible proportion by weight of the obscuring cloud. 

3. It is just these particles, however, which will be most influenced by 

the pressure of the radiation of the stars. Calculations from more ac- 
curate data confirm Schwarzchild’s conclusion that for a particle of the 
optimum size and the density of water, the repulsive force of the sun’s 
radiation is about ten times the gravitational attraction, and show that 
for stars of the same brightness, but other spectral types, the radiation 
pressure will be about two-thirds as great for Class M and increase for 
the whiter stars, till for Class B it is fully ten times as great as for solar 
stars. 
’ Dwarf stars will hardly repel dust at all, but giant stars, and especially 
the very luminous one of Class B, will repel it very powerfully. Only 
the coarser particles can come near such a star—the finer ones being driven 
away. ‘This selective removal, from the vicinity of bright stars, of the 
particles which are most efficient in cloud formation, may explain the 
fact that the luminous portions of these dark nebulae, though centered 
upon stars, do not brighten up in their immediate neighborhood as much 
as might have been anticipated. 

The finest dust must continue to be repelled by the stars, whatever 
their distance. It may congregate to some degree in interstellar regions, 


# 
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where the repulsive forces from stars on opposite sides are nearly equal, 
but it can be in no true equilibrium there, and must ultimately escape to 
an indefinitely great distance. 

4. Some force, however, operates to hold these dark clouds together, 
for their outlines are often sharp. This is probably the gravitational 
attraction of the cloud itself. 

Taking a spherical cloud as an example we find that, if its mass is M 
times that of the sun, and its radius R parsecs, the velocity of escape at 
the surface is 0.092M’*R-“km./sec. ‘The internal velocity of the neb- 
ular material is known only in the case of the Orion nebula, where the 
luminous gas shows irregular variations in radial velocity from point to 
point, amounting to about 5 km./sec.* on each side of the mean.- 

For a nebula 1 parsec in diameter (which may be taken as a rough rep- 
resentation of the small black, almost round spot about 15’ in diameter, 
discovered by Barnard® in Ophiuchus) the mass must be 60 times that 
of the sun, if the escape velocity is to be 1 km./sec. 

If all this matter were in the form of particles of rock of the optimum 
size, the extinction for light passing centrally through the cloud would 
be 2000 magnitudes. An extinction of 10 magnitudes (quite sufficient for 
opacity) would be produced if the radius of the particles were 72 yu. 

Though these numerical values are largely conjectural, it appears prob- 
able that the aggregate mass contained in one of these great obscuring 
clouds must be very considerable—probably sufficient to form hundreds 
of stars—and that a sensible fraction of the whole mass must be in the form 
of dust less than 0.1 mm. in diameter. 

It can easily be shown that any dust cloud which is impervious to light 
must also be impervious to particles such as those of which it is composed 
(and to free-moving electrons as well) in the sense that such a particle 
could not traverse the cloud without a practical certainty of collision. 
These collisions may account for the existence of dust within the clouds, 
even if it was not a primitive constituent. 

5. The transition from these dark nebulae to luminous nebulae in the 
vicinity of the stars appears to occur in two ways. ‘The first is by simple 
reflection of the light of the stars—which appears to occur in the nebulosity 
surrounding the Pleiades, the star p Ophiuchi, and probably in many other 
cases. ‘The second is by the excitation of gaseous emission, as in the Great 
Nebula of Orion, which is connected with one of the greatest known re- 
gions of obscuration and itself shows signs that obscuring masses lie in 
front of it. 

Both theoretical considerations, as suggested by the writer® and the 
facts of observation collected by Hubble,’ indicate that the luminosity 
of gaseous nebulae is probably due to excitation of the individual atoms 
by radiations of some sort (ethereal or corpuscular) emanating from neigh- 
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boring stars of very high temperature. In the Orion nebula the stars 
of the Trapezium (@ Orionis) appear to be the source of excitation. 

There is no reason to believe that the luminous gas forms the whole, 
or even any large part, of the matter present withm the region—only 
that it is selectively sensitive to the incident excitation, and therefore 
gives out most of the light, just as the gases (carbon compounds and nitro- 
gen) do in the coma and tail of a comet. 

If the turbulent motions of the various parts of this nebula are of the 
same order of magnitude in the other two coérdinates as in the radial 
direction, they must correspond to an average proper motion of 1.5 as- 
tronomical units per year, or about 0”.8 per century (with Kapteyn’s 
parallax of 0”.0055). In a million years this would carry a nebulous wisp 
through 2°, which is more than the whole extent of the nebula. 

It appears probable, therefore, that the aspect of the Orion nebula 
was entirely different a million years ago from what it is now, as regards 
its details. There is no reason, however, to suppose that the nebula 
was not there. We may rather imagine that wisps and clouds of dust, 
carrying gas with them, are slowly drifting about. Some of them pass 
through the field of excitation due to the radiations from the Trapezium 
stars, and, when in this field, the gas is set shining—faintly near its out- 
skirts, and without excitation of the nebular lines; more strongly, and 
with the nebular lines, near the middle. 

According te unpublished investigations by Hubble, it appears probable 
that the absorbing clouds in Orion, not far from the nebula, weaken the 
light of stars behind them by at least ten magnitudes. The exciting 
radiations probably penetrate to a relatively small depth into the mass 
and, even if they went deeper, little of the excited light could get out again. 
The Orion nebula, on this hypothesis, may be regarded almost as a super- 
ficial fluorescence of the gaseous portion of this vast dark cloud, in the 
limited region where it is stimulated by the influence of the exciting stars. 

1 Barnard, E. E., Astrophys. J., Chicago, 49, 1919 (1-23). 

2 Seares, F. H., and E. P. Hubble, Jbid., 52, 1920 (8-22); Mt. Wilson Contr., No. 187. 

3 Schwarzschild, K., Sitzungsberichte der K. B. Akad. der Wiss., Math. Phys. K1., 
Miinchen, 31, 1901 (293-338); Proudman, Monthly Not., R. A. S., London, 73, 1913 
(535-539). 

4 Barnard, E. E., Astrophys. J., Chicago, 38, 1913 (496-501). 

5 Publications of the Lick Observatory, Berkeley, Cal., 13, 1918 (98). 

6 Russell, H. N., The Observatory, London, 44, 1921 (72). 

7 Hubble, E. P., “Annual Report of the Mount Wilson Observatory,” 1921; 
Year Book of the Carnegie Institution of Washington, 1921. 








